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ABSTRACT 

A  new  method  for  predicting  the  off  design  performance  of  a 
multistage  compressor  has  been  developed.      This  method  replaces 
the  usual  routine  numerical  integration  through  the  compressor  by 
approximate  analytical  integration.      The  overall  compressor  perform- 
ance of  the  compressor  is  found  as  a  perturbation  of  the  design 
performance  in  terms  of  linear  and  quadratic  expansions  in  two 
parameters,   the  deviation  from  the  design  flow  rate  and  the  deviation 
from  the  design  speed.      The  procedure  for  this  preliminary  examina- 
tion has  been  simplified  by  assuming  that  all  the  single  stage 
performance  characteristics  are  identical  and  that  the  difference 
between  static  and  stagnation  conditions  are  negligible. 

The  results  of  this  method  are  quadratic  equations, in  terms  of 
the  two  perturbation  parameter s, for  the  temperature  ratio,   pressure 
ratio  and  adiabatic  efficiency  of  the  overall  compressor.      The 
amount  of  numerical  work  in  finding  the  coefficients  of  these  quadratic 
expansions  is  not  negligible,   however  the  method  appears  promising. 
Reasonably  accurate  results  are  expected  for  speeds  varying  as  much 
as  ten  per  cent  from  design  and  flow  rates  varying  as  much  as  four 
per  cent  from  the  design  flow. 
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I.     INTRODUCTION 

The  usual  technique  for  calculation  of  multistage  compressor 
performance  is  numerical.     Starting  at  the  entrance  with  a  prescribed 
flow  rate  and  rotor  speed  the  pressures,   temperatures  and  flow 
coefficients  are  determined  behind  each  stage  in  succession  from  the 
known  performance  characteristics  of  the  individual  stages.     Although 
this  method  can  be  carried  out  rather  rapidly  on  machines,    it  is 
desirable  to  have  available  a  method  that  is  more  suitable  for  explor- 
atory investigations  of  the  influence  of  stage  characteristics. 

A  new  method  is  proposed  which  replaces  the  routine  numerical 
integration  through  the  compressor  by  approximate  analytical  integra- 
tion.     The  overall  performance  of  the  compressor  can  be  found  as  a 
perturbation  of  the  design  performance  in  terms  of  linear  and  quadratic 
expansions  in  two  parameters,   the  deviation  from  the  design  flow  rate 
and  the  deviation  from  the  design  speed.      The  amount  of  numerical 
work  in  finding  the  coefficients  in  the  expansion  is  not  negligible, 
nevertheless  the  method  appears  promising. 

The  investigation  should  be  considered  as  preliminary.     Cer- 
tain approximations,    such  as  the  assumptions  that  all  single  stage 
performance  characteristics  are  identical  and  that  the  difference  be- 
tween stagnation  and  static  conditions  are  negligible,    are  not  essential 
but  simplify  the  procedure  for  the  initial  examination. 
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II .     DEVELOPMENT  OF  EQUATIONS 

A.  Work  Coefficient 

The  performance  of  multistage  compressors  is  generally 
presented  as  plots  of  temperature  and  pressure  ratios  against  dimen- 
sionless  weight  flow  for  given  dimensionless  compressor  speed.      The 
stage  performance  is  most  conveniently  represented  in  the  dimension- 
less  form  as  a  work  coefficient    (Y)>   pressure  coefficient    (Y)    and 
adiabatic  efficiency   (r\)    against  a  flow  coefficient    (4>).      These 
dimensionless  quantities  are  defined  as: 


Y 


A  h.  cfAT 

uvz   "  uvz 


r-  f& 


* 


A 

0 


U  pUA 

Y 


Y 


Ah       is  the  increase  in  stagnation  enthalpy  per  stage,     c       the  specific 
heat  of  the  gas  at  constant  pressures      AT    the  increase  in  stagnation 
temperature  per  stage,     Ap    the  increase  in  stagnation  pressure  per 
stage,   and    p    the  gas  density.      The  compressor  rotor  reference 
speed  is  represented  by    U    and  the  mean  axial  gas  velocity  by    V 
m    is  the  mass  flow  of  gas  and    A    the  cross  sectional  area  of  the 
flow  annulus. 
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The  multistage  compressor  temperature  ratio,   pressure  ratio, 
and     adiabatic  efficiency  will  be  determined  from  the  single  stage 
dimensionless  characteristics. 

In  general  it  can  be  said  that  the  work  coefficient    V     is  a 
function  of  the  flow  coefficient: 

Y    =   Fl-H 

The  temperature  rise     (AT)    in  a  single  stage  can  thus  be  written  as  a 
function  of  the  flow  coefficient: 

A  multistage  compressor  will  be  assumed  to  be  made  up  of 
stages  which  all  have  the  same  characteristic.     Such  a  compressor 
can  be  considered  to  be  made  up  of  an  infinite  number  of  infinitesimal 
stages  such  that  the  temperature  can  be  considered  to  be  a  continuous 
variable  along  the  axis  of  the  compressor.      Let  the  continuous  axial 
variable  be    £,     replacing  the  sequence    £       =    nt'   wnere    n    ^s  the 
stage  number  from  the  entrance  and    N    is  the  total  number  of  stages 
of  the  compressor     (0  ^  n  5.  N).      Using  the  variable    £    the  tempera- 
ture rise  can  be  written  in  the  differential  form: 

4r     -    —  t(>) 

The  compressor  then  has  an  overall  "length"  of  unity,    and  the 
function     T  (§)    becomes: 
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The  adiabatic  compressor  efficiency  is  normally  defined  as: 

=      T;    -     T<«> 

1        t  -  r<°> 

where    T    is  the  actual  temperature  and    T,     is  the  isentropic  final 

temperature  for  a  given  pressure  rise.     Introducing  the  isentropic 

Ti  r       f-J* 

relationship      „ ,    i  =  "M-t-j  this  becomes: 

r      T(o)       [p(oJ  ^ 


T/TG>)     —    I 

where     "/    .      and    T/T(o)    are  the  total  pressure  and  temperature 
p(o)  r  r 

ratios,      p(o)    and    T(o)    are  the  initial  pressure  and  temperature  and 

0     is  the  ratio  of  specific  heats  of  the  gas.     In  a  single  stage  the 
pressure  and  temperature  increase  can  be  represented  by    Ap    and 
AT     such  that    p    =    p(o)  +  Ap    and    T  =    T(o)  +  AT.      The  stage 
efficiency  thus  becomes: 

(i  +  *V^  -  I 


n   - 


AT 

Since    Ap    is  small  compared  to    p(o)    the  pressure  term  can  be 
expanded  to  give: 


¥(%)^m#p4+--- 


n    = 

AT 


Therefore  in  the  limiting  case,   with    Ap    very  small  (as  in  the  single 
stage): 


Iji     AP/pw 
A  T/T(o) 


or 


pw  (  1  ir-|  /    TCo)  (D 

This  will  also  be  the  form  of  the  limiting  case  as     Ap — ►  o     so  that  the 
differential  form  can  be  written  as: 

_L     Jfi.       _       77  /_XJ    _L    Jl_  (la) 

p  7T  u-i/  t    jc 

The  adiabatic  efficiency  (r|)  as  defined  here  can  be  shown  to 
be  equivalent  to  the  definition  r\  =  /  "ty"  ^n  the  limiting  case.  By 
the  definitions  of  the  work  and  pressure  coefficients: 

V  W7z)  AP 


Y  cpAT  o  Co  AT 

But  for     Ap     small  equation  (1)  applies  and: 

y  e  cp 

Substituting  the  thermodynamic  relations  for  a  perfect  gas: 

p  =  R  _R V-l 


CP 


gives  _  * 
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B.  Differential  Equations  for  Flow  Coefficient 

The  flow  coefficient  has  been  defined  as: 

9     =       pUA 
Introducing  the  perfect  gas  relation, 

the  flow  coefficient  becomes:  _ 

,  rh  K   / 

Taking  the  logarithm  and  differentiating  with  respect  to    £    this 
becomes: 

4'    _      r/        p'  A' 

<f>  T  >  A 

where  the     (     '  )    denotes  differentiation  with  respect  to  the  variable  £. 

p' 

Substituting  for     —    from  equation  (la)  this  reduces  to: 

.*'     -    /■-  n  -*-}   -XI  _  _Al  (2) 


f 


4>  L  *-'J      T  A 

This  is  an  approximate  relation  since  the  differences  between  static 
and  total  temperature  and  pressure  have  been  neglected  in  the 
development. 

At  the  compressor  design  point  the  flow  coefficients    efficiency 
and  temperature  rise  are  assumed  to  be  constants  throughout    £,    and 
equal  to  the  single  stage  parameters  at  the  design  point.      Let  the  sub- 
script   (    0)    denote  the  design  condition.      Therefore  at  design 
conditions: 


<f>     =    cj)       =    constant 


r\     =    n       =    constant 


*     =   T'=    — 1=    constant 


(3) 


Vd4  A       IO       c, 

At  the  design  point  equation  (3)  can  easily  be  integrated  to  give  the 
temperature  at  any  point  in  the  compressor  as: 

T.   —  r*  €     +    Oo  (4) 

where       @„    is  the  inlet  temperature. 

Dividing  equation  (3)  by  (4)  gives: 


(5) 


U.1  T. 

T. 

a 

€  ♦£ 


where    3       is  defined  as: 

1  U    r. 

Equation  (2)  can  now  be  solved  to  determine  the  area  variation 
of  the  compressor: 


1_    = 
A 


ftH 


nJ-^-r)  -4-£l  -  -A- 


Letting  the  term    r|(  .  _    )  -  1     =    f    and  noting  that    c|>     =0     (since    cj) 
is  a  constant)  and  substituting  in  equation  (5): 

A  =       _  ^  (6) 

A  ?    +   (9. 
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This  area  distribution  is  fixed  by  the  design  geometry.     Substituting 
(6)  into  equation  (2)  the  general  differential  equation  becomes: 


£ 


4> 


-  f 


r 


T 


+■ 


f. 


C  +  (3. 


(7) 


To  eliminate  the  temperature    (T)     equation  (7)  is  solved  for    T    and 
substituted  in  the  expression  for  T" : 


r  = 


ft 


r  = 


u 


fc-r 


r     — 


J 


u*    jc 


f-^rJ-i- 


€+(*. 


_  _4 


i. 


f  ^ 


c,T' 


u: 


J4 


•f  r 

- 

1       f- 

4>     J 

Performing  the  differentiations  and  substituting  the  relations: 

<Jr  Jr     d4>  dr 


df 


d<(>     dS 
df      d4> 


d<£ 
df 


* 


d£  d<f       d£  d<f 

there  results  the  general  differential  equation  in    <{>: 


fr<j>" 


fs_  ffrif 

+  0.  LA    d(j> 


4*0 


4-  f 


d<f> 


+  ^r 


']f 


(8) 


-  IMI  *'«♦"*•#♦  J#3£t 


O 
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C.  Solution  for  the  Flow  Coefficient 

A  series  expansion  of  each  of  the  terms  involved  in  the 
differential  equation  has  been  used  to  find  a  solution  for  the  flow 
coefficient.     From  the  single  stage  performance  it  is  seen  that  both 
the  efficiency  function      f    and  -work  coefficient      V     are  functions  of 
the  variable    <\>: 

f  -  f  (♦)  r  -  r(<f>) 

These  functions  can  thus  be  expressed  as  a  Taylor  series  expansion 
about  the  design  point.     The  stage  efficiency  at  design  is  assumed  to 
have  its  maximum  value,    hence  the  first  derivative  of  the  efficiency 
function  becomes  zero     ( ,     ■)      =  0.      The  Taylor  series  expansions 
thus  become: 


™-  T-  +  (#<♦-*!+  zfeMn  *{#)(♦-♦£- 

The  flow  coefficient  will  be  represented  by  a  series  of  the 
form: 


(9) 


:io) 


The  factor     £    represents  a  convenient  small  parameter  to  be 
replaced  by  unity  in  the  final  form.     Substitution  of  the  series  (9)  and 
(10)  into  equation  (8)  results  in  a  series  equation  in  powers  of    €.  . 
Equating  terms  in  like  powers  of      €.     gives  a  series  of  linear 
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differential  equations  in  the  flow  coefficients    4>     (n  =  1 ,    2  .  .  .  .  n). 
Only  the  first  two  of  these  coefficients  will  be  retained.      The  linear 
differential  equations  for  these  coefficients  are: 

4>;  +■  {*£♦.  ^i}^j:-  o  (id 


f+{ti*  ^f  d*. 


♦  tl  4>, 


f.  J  ♦. 


+K^- #-*(?♦-&+■-* 


Ml 

Z+p.  (12) 


Z    T*    £+A_    r« 


S+fr 


/ 


Note  that  the    (        )    denotes  differentiation  with  respect  to  the  varible 
£    when  used  with  the  flow  function    cj>    and  denotes  differentiation  with 
respect  to  the  variable    <j)    when  used  with    f    or    T*  .     The  constants 
f    ,     1"'   ,     T     ,   f    ,    f  and    <j>       are  all  known  from  the  single  stage 

characteristic. 

To  simplify  the  solution  of  equations  (11)  and  (12)  the  following 
notations  have  been  introduced: 

r.    T* 
oc    =     *.  +  £ 

*  r.     +  X 
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The  coefficients  of  the  left  side  of  equation  (12)  become: 

B=  r^{5-  +  ^^'}   ;     c-^'(?.-f.)-U*t;J)-ta 

Using  this  notation  equations  (11)  and  (12)  become: 

+"  +    *?£*'-   °  (lla) 

<b"  +        -^<f>'  =       B4>V     C      A£    +  £_&!_  (12a) 


Multiplying  equation  (11a)  by  the  integrating  factor     (£  +  (3   ) 


<1€ 

Integrating: 


(te+zO*  ♦,'  }  -  o 


(£,  ¥  (5a)       (t>(        =    constant 

The  constant  can  be  found  by  introducing  the  boundary  condition 
at    %     -    0;    (J),     =    <j>,    (o)    and  therefore: 

Constant    =       /3a       (J^  (o) 


Therefore: 


4,'=  (r^arf  <*)''<0)  <Ub> 
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Integrating  between  the  limits    0    to    £: 


(o) 


In  terms  of     £ft     this  becomes: 


(j)(      =     <j)((o) 


0. 


(5. 


?.+"!    IVC  +  p- 


*.+' 


4>;<o) 


(13) 


Equation  12a  is  solved  in  a  similar  manner.     The  integrating  factor 
(£  +  (3    )        is  again  applied  to  the  left  side  of  the  equation  and  the 
solutions  for    §,     and    $.     form  (lib)  and  (13)  substituted  into  the 
right  side.     On  integrating  once  the  equation  is  of  the  form: 


fc'-  Eft- P-) 


t/^.  +  3) 


F(*-»-j3.)      +G(S+p.)   +H(|+p.)         Inl±£ 


The  coefficients  E ,    F,   G,    and    H    are  a  collection  of  all  the  constants 
in  the  equation  not  involved  in  the  integration.      On  integrating  the 
second  term  of  the  flow  coefficient  becomes: 


*>  -  *«  -  £iM&-.)  '-'J- F^  feW  -  ■} 


(14) 


The  coefficients  of  equation  (14)  are  defined  as  follows: 
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2iX.+l) 


E  =  -{B-£r  +  c^}f7T     4>« 


'  I        t.  +  t)A.  +  t     r>  (?.w)»  r*     ti      T, 


(15) 


All  terms  containing  the  variable    £    will  be  designated  by  lower  cas. 

b's  with  subscripts: 

f       A     \2<*  +  '        "I 

ki  =  sSil(rw  "  'J 

"-=  -frfiAT-  -] 

c+olU+p./       'J 


2(X.+01U  +  (W  <16) 


bs-     In 


fr 


€*?« 
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Introducing  the  coefficients  (15)  and  (16)  into  equations  (13) 
and  (14)  the  first  and  second  flow  coefficient  terms  become: 

(j)(   -     ^(o)     _      b3   <j)/(0)  (13a) 

(j)2_      (fc  (0)    4.  Q   (K(o)  4-C^0»)  <j>(o)4-d  (j)/o)      4-  e  (j)z<°)  (14a) 


where  the  coefficients  of  (14a)  are  defined  as: 

Since  a  large  number  of  coefficients  will  be  used  in  the 
following  equations,    they  will  be  represented  by  lower  case  letters 

with  number  subscripts.      The  letter  will  be  used  to  indicate  coeffi- 

2  /  2 

cients  of  a  given  term  (a    for    <)>,  (0)5  c    for    <j>,  (o)    §,  (o);    d    for    4>,(o) 

and  e  for  either  <(>,  (o)  or  <j>?(o)).  The  number  subscript  refers  to 
a  given  equation.  The  lower  case  b  refers  to  terms  containing  the 
variable    £. 

D.  Flow  Coefficient  Initial  Conditions 

Equations  13a  and  14a  can  be  used  to  determine  the  flow 
coefficient  in  terms  of  the  initial  conditions  at    £     =    0;   (4>,  (o),   <f>,  (o), 
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cj)    (o)    and    (j)_(o)).      The  first  stage  of  the  compressor  will  be  assumed 
to  operate  with  the  same  characteristic  as  the  single  stage,    therefore 
at    £     =    0    A(o)    will  be  known.      Let    §{o)     =    <(>     +  e<t>,(o)    with  (j>    (o)  =  0 
for    n    >  1,   thus    ^ifa)    is  prescribed  and    (j>?(o)     =    0.     For  the  initial 
condition  at      £     =    0    equation  (7)  becomes: 

A'w_  (f(o,  _Tj£>     +     _L_)    ((,(0)  (17) 


By  definition: 


'dTV 

> 

r(o) 

rw  = 

o 

u* 

r(o) 

u 

*  r„ 

(18) 


T(o)  u   ^(o)       u*  rc     /  r(<0 


T^o)  cp  o  cr<e>    V   T. 


The  term     —  is  similar  to  the  term  defined  as  yC.  ,   therefore  let: 

Substituting  from  the  expansions  (9)  and  (10)  gives: 

T(°)  —     ^-  -f-€T.'4>,<0)-h  e*  -?H  cbcof 

{to)  -     f.  +  e'4-4,'(0)Z  (20) 
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Substituting  (18)  and  the  expansions  (19)  and  (20)  into  equations  (17) 
results  in  a  series  equation  in  powers  of    £    .     Equating  terms  of  like 
powers  of     £     gives  the  following  expressions  for  the  flow  coefficient 
derivatives  at  the  initial  conditions: 


f>=  -ir*   -  -fr*.'" 


(21) 
z 


**-  {*&-)*•*•-  J#-'-  iK +  *}♦> 


E.  Flow  Coefficient  -  Final  Form 

Substituting  the  initial  conditions  (21)  into  the  expressions  for 
the  flow  coefficient     terms  (13a  and  14a)  results  in  a  quadratic 
expression  for  the  flow  coefficient  in  terms  of  the  velocity  (or  speed) 
perturbation       (3,  ,   and  the  flow  coefficient  perturbation       $,  (o)    of 
the  form: 

4>(4)   =    ((>. +  eA1p,  +  eM,^,VA«^+€  A,  (()/•)+ eM,,^)1 

The  upper  case  letters  will  be  used  to  represent  coefficients  of  a 
given  equation  and  the  subscripts  refer  to  the  terms  in  the  equation. 
The  subscripts  1  and  2  will  be  used  to  represent  coefficients  of  the 
terms    §■.  (o)    and    p.     respectively. 

Small  deviations  from  the  design  point  are  accomplished  by 
letting  £  4>.  (o)  and  £  (3,  be  small.  This  can  be  accomplished  by 
letting     €•  =     1     and  using  small  perturbation  of    <}>,(o)    and    p\. 
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Letting      6=1     the  flow  coefficient  becomes: 


7. 


The  coefficients  of  equation  (22)  are  expressed  in  terms  of  the  initial 
design  conditions  (•"•)'*•  ,'•  )T»jTo   >  <P»  )  and  the  location  along  the  length 
of  the  compressor    £.     These  coefficients  are  most  easily  expressed 
in  terms  of  the  coefficients  of  all  the  previous  equations  involved. 
Express  the  initial  conditions  (21)  as: 

cj>>)=    B,  <)>,(o)     +BJ, 

(21a) 


■where 

Substituting  into  equations  (13a)  and  (14a)  the  coefficients  of  equation 
(22)  become: 

Az  =  —Btb, 

/\it  -     a  B2*     -  b3B2i 

Alt-  Za3,B£  4-cBz     -baBlt 

/\(  -  I       -    b3B, 

y^     =      Q  B,   4-  cB,    4-  6       —  b3Bu 
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F.  Temperature  and  Pressure  Ratios 

The  differential  equation  for  the  temperature  rise  through  the 
compressor  is  integrated  to  determine  the  temperature  ratios  as 


follows:  (JT  U  y  (_<b) 

Tt5)-    T(o)    -      «yl  f   T(0   d€ 
Substituting  the  expansion  of     V  (§)    and  dividing  by      T(o)     =    0>   the 
temperature  ratio  is: 


™>      .     ,+ 


T(o) 


^ [%r.  +. « r/jVs  +■  c't.^  j?  +  ^*JJ 


U*  T.         i 

Let      £  =    1     and  substitute    —    =    -3-    the  temperature  ratio  becomes: 

'      tS  -  '  *■  t{«  +  # + *JJ?  +  « i*+: J§)     (23) 


The  flow  coefficient  is  defined  as: 
The  flow  coefficient  ratio  at  any  station  in  the  compressor  becomes: 

4>(g)         m R T(g)      phU^A")  _     rcc)  pc^i    /U») 

4>(o)  rnRTW        p<*>U<«>A<*)  *n°)    pc?)     A(0 

since    m    and    U    are  constant  through  the  compressor.      The  pressure 
ratio  is  thus  given  by: 

pit)  T(5)    4><°)        A  (•) 


co)  T(o)     4>(S)       A  <€) 


(24) 
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The  temperature  ratio  is  given  by  equation  (23)  and  the  flow  coeffi- 
cients by  equations  (20)  and  (22).      The  area  ratio  is  fixed  by  the 
geometry  at  the  design  condition  given  by  equation  (6)  as: 

i       dA  ft. 


A       dS  VS  +  P. 

Integrating  gives:  f  e 


I 


t 


A  'H  +  f 


,h  a^t  -    f- ,h  t^; 


AW 


A  (5) 


f. 


The  integrals  in  equation  (2  3)  are  found  by  direct  integration 

2 

of    (J),  ,   (J>,     and    <j>?     as  given  by  equations  (13)  and  (14).      The  results 

are  again  expressed  in  terms  of  the  initial  conditions  and  coefficients 
involving  the  design  conditions  and  axial  position  (£)    through  the 
compressor.     Thus: 

A  ,, 

>  2  2. 
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In  terms  of  the  coefficients  of  equation  (12)  and  the    £    variable 
coefficients  (16)  the  coefficients  of  these  equations  are  given  by: 


e—  e, 


Substituting  these  integrals  into  equation  (2  3)  and  introducing 
the  initial  conditions  as  given  by  equations  (21a)  the  temperature 
ratio  also  reduces  to  a  quadratic  expression  in  the  velocity  perturba^ 
tion    p,     and  the  initial  flow  coefficient  perturbation    <j)1  (o).      This 
equation  is  of  the  form: 
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t^t"  ~k[i^^) ^  ^tf  +  J)«^*»(d)+:d' ^ ^^"^ 


.     (23a) 


The  coefficients  of  equation  (23a)  are  given  by: 


2± 
-a.      r.  ~z 


U-  ^Bte,    - 


Dl  -  4  (  B,eI+  ?) 

D,  -  ^(zq.Bft  +  B^.^e)*  -£(2«£B,B,+-  B2ca)_  -£_ 

n-   -£l  (B,a  a3  +  3,  cs  +  d  3  +B.e.)  +  j£  (aaB,V  B,  c^  ?) 


The  maximum  value  of  the  temperature  ratio  at  a  given  speed 
({3,  =  constant)  is  found  by  setting  the  derivative  of  the  temperature 
ratio  with  respect  to    <j>,  (o)     equal  to  zero: 

_i_  {JB£,}_   -L_(d  S.-hD,    4-21>,<j>4==0 

The  value  of  the  flow  coefficient  for  the  maximum  temperature  ratio 

is  thus: 

D,x  (3,  +■  D, 


(f)(0)  =_      

•  I        max 


zvtl 


The  second  derivative  corresponds  to  the  curvature  of  the  temperature 
ratio  curve  and  is: 
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J z       {  T(j)\  _  ZP„ 


and  is  independent  of    p\.      The  temperature  ratio  curves  can  thus  be 
completely  determined  by  the  coefficients  of  equation  (23a). 

A  similar  parabolic  curve  can  be  found  for  the  pressure  ratio 
by  substituting  equations  (20),    (22),     (2  3a)  and  (2  5)  into  equation  (24), 
and  dropping  all  the  terms  of    (3,     and    <f>,  (o)    higher  than  the  second 
order.      The  same  procedure  can  be  used  to  express  the  adiabatic 
efficiency    r\    by  substituting  the  equations  of  pressure  and  tempera- 
ture ratio  in  the  equation: 


V    - 


T(g) 


The  compressor  efficiency    (n    )    is  determined  by  setting    £     =    1. 
The  second  order  equations  for  the  pressure  ratio  and  efficiency  thus 
reduce  to  the  following  form; 

-pfig  =  (11&.)  '*  (l  4-  E£  ft  +  Ei2  (3,2  -.-  Eaft4>.M-  £,  ♦,*+  E„+>j}    (25) 


>?  -  v.  (_  i  +■  c*  p.  -•-  c«  ?!+c.2ftt/«i+  c,  4,/«h-  a  <f>,">J 


(26) 


The  coefficients  of  these  equations  are  most  easily  expressed  in  terms 
of  the  coefficients  of  the  previous  equations  and  are: 
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r  T>%  _      _Ai  £  I  -A,     +  J>, 


^     Ifcl         4>.         4>.Ce+p.)        £  +  p. 

r         f  A'f  4 4lL    +    -J—(d    +     D'-^'P') 


1         K-l 


r K_    ]M/  F  _FD £*-)  +  D2  -  D 

Or    K-l         If     tfc«       L*U*         2*J    +^         ^ 


The  coefficient    C,     can  be  assumed  to  be  zero  since  the 
maximum  efficiency  at  design  speed  is  known  to  occur  at  the  single 
stage  design  flow  rate.     Setting  the  derivative  of  the  efficiency  with 
respect  to  the  initial  flow  coefficient  equal  to  zero: 

the  maximum  efficiency  is  found  to  occur  at  an  initial  flow  coefficient 
perturbation  of: 
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*•-  -  {-^fe-} 


Since  the  maximum  adiabatic  efficiency  at  design  speed    Q3     =    0)    must 
occur  very  close  to  the  design  flow,     C,      should  be  very  small. 

The  speed  function    (p    )    for  maximum  efficiency  is  found  by 
substituting  the  value  of  the  flow  coefficient  perturbation    (cj).  (o))    for 
maximum  efficiency  into  equation    (2.6)    which  becomes: 


^_^|+(Cl-^f)(3.   +(C-     -^)(3 


Setting  the  derivative  with  respect  to    (3,     equal  to  zero  the  maximum 
efficiency  is  found  to  occur  at  a  speed  function  of: 


p, 


C,  C,t     -    Z  C„  Cz 
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III.     APPLICATION 

The  performance  of  a  ten  stage  compressor  has  been  computed 
to  illustrate  the  use  of  the  method.      The  single  stage  characteristics 
for  the  research  compressor  as  described  in  reference  1  have  been 
used  in  this  analysis.      The  single  stage  characteristics  for  this 
compressor  are  shown  in  Figure  1  in  terms  of  the  parameters      T 
and    f.      The  values  of  the  design  point  constants  for  this  example  are; 

r    =  1.95,    r'   =  -7.25,    r"  =-37.5,   f    =2.01,   f"  =  -  1  is 

00  000 

and    <})       =0.  44.     The  reference  rotor  speed  for  this  example  is 

U       =    1000  fps.      The  gas  is  air  at  standard  temperature  and  pressure 

with     Q        =    520°R,     C       =    .  24    Btu/lb°R     and      tf  =  1.4.      The  value 
o  p 

of    3       and     z;       then  become;     (3       =     1.603    and     t,        -  -  3.288.      The 
ro  o  ro  o 

broken  line  curves  in  Figure  1  are  the  approximate  curves  represent- 
ing the  quadratic  expansions  found  by  dropping  all  terms  higher  than 
the  second  order  in  equations  (9).      This  approximation  is  seen  to  be 
reasonably  accurate  for  flow  coefficients  from  0.  40  to  0.  50. 

The  overall  compressor  performance  is  found  by  letting 
N     =    10    and    £     =    1 .  0    in  equations  (23a),    (25)  and  (26).      The  speed 
variable  is  computed  from  equation  (19)  which  can  be  simplified  to 
give: 

The  resulting  temperature  ratio  curves  for  the  ten  stage  compressor 
are  shown  in  Figure  2  as  a  function  of  the  inlet  flow  coefficient.  The 
pressure  ratio  and  efficiency  curves  are  shown  in  Figure  3.      The 
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maximum  efficiency  is  found  to  occur  at    (3,     =    0.  363  which  corre- 
sponds toU/U       =    0.  903. 

The  operating  point  of  each  stage  is  found  by  varying    £     in  ten 
equal  increments  and  solving  equation  (22)  for  the  flow  coefficient. 
The  operating  points  at  design  speed  for  the  first,    second,    fifth  and 
last  stage  of  the  ten  stage  compressor  are  shown  in  Figure  4.      The 
stage  flow  coefficient    (<j>)    is  plotted  as  a  function  of  the  compressor 
inlet  flow  coefficient    (cj>  (o)).      The  stage  efficiency    (r\) ,    work  coeffi- 
cient    ("40    and  pressure  coefficient    (ly)    are  shown  as  functions  of 
the  stage  flow  coefficient    (§).      These  curves  indicate  the  range  in 
which  the  quadratic  approximation  is  reasonably  accurate.      The  quad- 
ratic expansion  of       i      and    f    are  shown  in  Figure  1  to  be  reasonably 
accurate  for  stage  flow  coefficients  between  0.  40  and  0.  50.     The  tenth 
stage  is  seen  to  be  outside  this  range  at  design  speed  for  compressor 
inlet  flow  coefficients  greater  than  0.  455  and  less  than  0.  400. 

The  last  stage  flow  coefficient  will  have  the  greatest  variation 
from  the  compressor  inlet  flow  coefficient  for  off  design  speed 
conditions.      The  operating  points  for  the  tenth  stage  at  varying  speeds 
are  shown  in  Figure  5. 

The  operating  range  of  the  tenth  stage  indicates  that  the  method 
should  be  reasonably  accurate  for  a  fairly  wide  range  of  compressor 
inlet  flow  coefficients  at  design  speed.     For  varying  speeds  reasonable 
accuracy  can  be  expected  at  speeds  varying  up  to  ten  per  cent  off 
design  with  a  similar  flow  coefficient  deviation. 
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CONCLUDING  REMARKS 

It  has  been  shown  that  the  off  design  characteristics  of  a  multi- 
stage compressor  can  be  predicted  by  a  method  using  approximate 
analytical  integration.      The  method  can  be  expected  to  give  reasonably 
accurate  results  at  design  speed  for  compressor  inlet  flow  coefficients 
varying  as  much  as  four  per  cent  from  the  design  condition.     At  off 
design  speeds  the  method  should  be  accurate  for  speeds  varying  up  to 
ten  per  cent  from  design  with  similar  flow  coefficient  variation  for  each 
speed.     A  comparison  with  actual  test  results  has  not  been  made  since 
none  could  be  found  for  a  compressor  with  all  stages  having  identical 
characteristics. 

This  method  should  be  extended  to  include  varying  stage 
characteristics.      This  could  be  done  by  assuming  a  given  single  stage 
characteristic  over  a  section  of  the  compressor  and  different  single 
stage  characteristics  over  other  sections.      The  approximate  integra- 
tions could  then  be  carried  out  over  the  various  sections  using  the 
appropriate  stage  constants  for  that  section. 
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Fig.   2         Ten  Stage  Compressor  -  Overall  Temperature  Ratio. 
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Fig.   3         Ten  Stage  Compressor  -  Overall  Pressure  Ratio 
and  Adiabatic  Efficiency. 
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